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Abstract 

For physical theories, the degree of arbitrariness of a system is of great importance, 
and is often closely linked to the concept of degree of freedom, and for most systems 
this number is far from obvious. In this paper we present an easy to apply algorithm for 
calculating the degree of arbitrariness for large classes of systems formulated in the language 
of moving frames, by mere manipulation of the indices of the differential invariants. We 
then give several examples illustrating our procedure, including a derivation of the degree 
of arbitrariness of solutions of the Einstein equation in arbitrary dimensions, which is vastly 
simpler than previous calculations, the degree of arbitrariness of gauge theories of arbitrary 
groups under Yang-Mills type equations, for which until now only specific cases have been 
calculated by rather messy means, and the degree of arbitrariness of relativistic rigid flow 
with or without additional constraints, which has not been derived before. Finally we give 
the proof of our algorithm. 
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1 Introduction 



1.1 The question of the degree of arbitrariness of a system 

The degree of arbitrariness of physical systems is defined as the number of arbitrary func- 
tions that one must specify in order to completely determine the state of the system. 
Often it is also closely linked to the concept of the degree of freedom of a system, which 
is heuristically described as the number of ways the constituents of the system can move 
independently: for example, the degree of freedom in the kinetic theory of gases, or of 
elementary particles. The concept of degree of arbitrariness is general and not restricted to 
microscopic systems nor systems constructed with classical particles: we also talk about the 
degree of arbitrariness of a solid body, of a fluid flow, of the electromagnetic field, of a gauge 
field, or even of the gravitational field. As physical theories are written in terms of equa- 
tions, the importance of the degree of arbitrariness, interpreted as the number of arbitrary 
functions entering the solution of the equations of the system, cannot be overemphasised. 

However, despite the great importance of this quantity, there is so far no general strategy 
for calculating it. In the classical theory of gases or of rigid bodies, calculation is done 
mostly "by intuition" , and in slightly more complicated circumstances, ad hoc arguments 
and calculations arc used to suggest a degree of arbitrariness (feasible since one is almost 
always restricted to dimension 3 + 1), which is then justified a posteriori by checking 
that the suggested number is physically sensible. In the theory of fields, the degree of 
arbitrariness, when it is mentioned at all, is almost always calculated "by definition": 
the geometrical space of the theory itself is constructed using a number of free functions. 
These strategies work, either because the system is fairly simple, or the construction of the 
system is straightforward. But even in the latter case there are often many pitfalls: for 
example, the metric of spacetime has 6 degrees of arbitrariness, but a Riemannian geometry 
is constructed using a bilinear form, which involves 10 free functions. The usual argument 
goes that we must take into the "diffeomorphism invariancc" of spacetime into account, 
which subtracts 4 from the counting. Of course, such arguments are a bit hand-waving, 
and whenever a geometrical setting is defined, instead by postulating something like the 
metric, by specifying a connection, such arguments do not apply. Related to Riemannian 
geometry, we also talk about the degree of arbitrariness of Einstein's gravity, which is, 
roughly speaking two times the degree of freedom of the graviton: remember that when 
calculating the degree of freedom of something, we are really restricting to a Cauchy surface 
for the equations of motion, and a function together with its time derivative, which are 
really two independent functions on the surface, counts as a single degree of freedom. 

In more complicated systems, especially in systems that are defined by using one or 
more constraints, the above-mentioned strategies fail miserably. Consequently, when faced 
with such a theory, instead of having a clear picture of what the physical variables are, 
we deal with a large number of dependent variables linked together by interrelated partial 
differential equations. One can still do calculations with such systems. However, due to 
the ignorance of the independent quantities in such systems, sooner or later calculations 
will have the tendency of degenerating into the random walk of some drunken sailor. 

In this paper, we propose an algorithm for calculating the degree of arbitrariness of a 
system by geometrical considerations. Our algorithm works by first formulating the phys- 
ical problems using the method of moving frames. As physical theories are often generic 
in nature: functions are not explicitly given, and are required to have their invariancc 
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properties: two systems are considered to be equivalent if a allowed transformation trans- 
forms one to the other, the use of the method of moving frames immediately grants us an 
incredible advantage to our algorithm: as the method of moving frames uses no coordinates 
whatsoever, the degree of arbitrariness we obtain is the true, physical degree of arbitrar- 
iness: there is no business of things like "diffeomorphism invariance of the coordinates", 
or "gauge invariance", etc., as we shall see later 0. It is hoped that the proposed method 
will promote the method of moving frames, which has hitherto seen various physical ap- 
plications in formulating the kinematics of physical systems, to be used also in studying 
the dynamics of physical systems, since it provides a way of studying equations of motions 
and constraints without specifying any coordinates. 

As physical theories are ultimately formulated using differential equations, any calcula- 
tion of degrees of arbitrariness amounts to the study of solutions of differential equations. 
Our algorithm is no different, and at its roots it relies on the Cauchy-Kowalewski existence 
theorem for systems of partial differential equations. However, that everything we do is 
reducible to the Cauchy-Kowalewski theorem does not mean that once we know this the- 
orem we know everything about calculating degrees of arbitrariness for physical theories. 
There are two problems: first, as we have mentioned, physical theories are often given in a 
generic form (which may even have the dimension of spacetime as a variable), and without 
the explicit functional forms of the variables, the Cauchy-Kowalewski theorem is not dir- 
ectly applicable; second, even if a theory is explicitly constructed, a direct application of 
the Cauchy-Kowalewski theorem would involve too much work. Our algorithm, on the 
other hand, makes use of an already non-trivial consequence of the Cauchy-Kowalewski 
theorem, the Cartan-Kahler existence theorem, to devise a simple procedure that can be 
applied without understanding any of the underlying mechanisms. If nothing else, sim- 
plicity should be the merit of our algorithm: one is required only to know how to do 
index manipulations and apply simple reasonings in combinatorics in order to apply our 
algorithm. 

1.2 The historical origin of our algorithm 

Our algorithm is based on the method of moving frames, the Cartan-Kahler theory of 
exterior differential systems, and to a lesser extent, the method of equivalence, all of which 
are due to the French mathematician Elie Cartan. Indeed, Cartan himself became briefly 
interested in the problem of calculating the degrees of arbitrariness of physical equations, 
and specifically those of the equations of general relativity and Einstein's (failed) unified 
field theory. With the equations explicitly given, Cartan of course knows how to calculate 
the degree of arbitrariness (see the historical note on page [T8|) , and he communicated the 
methods he used to Einstein in a series of letters. At first suspicious, Einstein finally 
understood Cartan's methods, and he wrote in his reply to Cartan [15]: 

Ich habe Ihr Manuskript gelesen und zwar mit Begeisterung. Nun ist mir alles 
klar. Sie sollten diese Theorie ausflihrlich publizieren, denn ich glaube, dass sie 
fundamental Bedeutung hat. 

^It should be noted that the method of moving frames is not the same as the tetrad method widely used 
in physics: the tetrad method depends on coordinates: it just derives a new set of coordinates from old ones. 
The method of moving frames, on the other hand, use differential forms, instead of coordinate functions, as the 
fundamental vehicle for doing calculations. 
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(I have read your manuscript — enthusiastically. Now everything is clear to 
me. You should publish this theory in detail; I believe it is of fundamental 
importance.) 

At Einstein's suggestion, Cartan did publish his methods with physicists in mind as his 
audience [7]. However, perhaps for fear that the physicists will not be too comfortable with 
his method of moving frames, which is still being actively developed by Cartan in private 
at the time, he chose to formulate the method so that it applies to first order partial 
differential equations, and due to the arbitrariness of the coordinates (and for modern 
applications, gauge "degrees of arbitrariness") involved, and the fact that his formulation 
requires knowledge of all the relations among the quantities in the differential equations, 
which are easy to be missed out unless one is using the method of moving frames, his method 
is difficult to apply, and even Einstein made various mistakes in applying the method [15] . 
After Einstein gave up his unified theory of fields, the Cartan-Einstein correspondences 
stopped, and Cartan never again attempted to communicate his methods to the physical 
community. It is very sad his published method did not spread very far among physicists, 
both due to the difficulty of the method and Einstein's later isolation from the mainstream 
physics community. 

Our algorithm can be considered a modern and much easier to apply reformulation 
of Cartan's original methods. The key improvements made to the original method is the 
consistent use of moving frames and the theory of differential invariants, both of which are 
also developed by Cartan. Our method also applies in cases where the mere form, instead 
of the explicit dimensions, etc., of a physical problem, are given. 

In Cartan's days it is still quite rare, and usually deemed unnecessary, to think in higher 
dimensions or arbitrary dimensions. In his last book [9], Cartan gave many examples of 
applications of his theory of exterior differential systems to obtain degrees of arbitrariness 
of many geometrical problems by direct application of the Cartan-Kahler theorem, some 
of them quite relevant for physics. But since most of these problems takes place on familiar 
dimensions (dimension 2 or 3), there is not too much motivation for generalisations or trying 
to find common features in the calculations. It is perhaps also for this reason that Cartan 
did not even attempt to formulate an easier method for calculating degrees of arbitrariness, 
and there is no doubt that he would have done much better than the present author did if 
he lives to see this day. 

1.3 Structure of this paper 

After a brief discussion of the necessary settings, we will immediately state our algorithm 
in section [2] without proof. Then in section [3l we give several example applications of our 
algorithms. Several examples chosen are such that the degree of arbitrariness is well-known, 
as a check of the correctness of our algorithm. Historical notes are included with these 
problems. The examples also include some new results, including the degree of arbitrariness 
of a general Yang-Mills type theory, and the degree of arbitrariness of relativistic rigid fiow 
under various conditions. For more involved applications of the algorithm to more complex 
problems, see the papers by the author [20l[2T]. Finally, in section [3] we give the proof of 
our algorithm. 

The rather unusual and not very logical ordering we chose for our presentation is due 
to the consideration that the algorithm itself is simple enough to understand and apply to 
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physical problems without too much additional mathematical baggages, whereas the proof 
involves the Cartan-Kahler theory on exterior differential systems, which few physicists 
know about. The mathematical constructions and theorems on which our algorithm is 
based on are reviewed in two appendices in this paper: we choose not to point simply to a 
book of reference since some of these books are rather dense and have to be read in a linear 
fashion. In addition, the modern approach to the theory is often highly formal, whereas 
we prefer a more down-to-earth approach suitable for calculations. As a consequence, our 
presentation of the Cartan-Kahler theory and of the equivalence method is closer to the 
original presentation of Cartan [5l[6l[9] than to what is usually done today. 

2 The algorithm 

We will first state our algorithm. Our description aims to be precise, and as a consequence 
the statements in the following may seem long-winded and complicated. This is only 
because we try to avoid "heuristics" in our statements, and as our examples will show later, 
once a problem has been formulated in the language of moving frames, the difficulty of the 
applications of our algorithm is only that of rather simple problems in the combinatorics 
of the indices. 

2.1 The moving frame 

As we have mentioned, first of all we must place the system we want to study into the 
language of moving frames. Assume this has been done: thus let cj^, n = 1,2, . . . , N be 
a co-frame for a suitable manifold P constructed from the geometrical problem. We can 
then proceed to write down the structural equations 

(2.1) dUfj, = Ifj,u\^i, A ux, 

where we have employed the Einstein summation convention. The quantities If_iy\, which 
are scalar functions on P, are called the fundamental invariants of the system. In addition, 
we may have a set of functions Ja subject to certain constraints. We further assume that 
we can find a subset of the 1-forms, denoted by Wj, i = 1, 2, . . . , n, such that all o;^, I^^x and 
Fa are labelled using indices i,j,k... instead of fi,u, p . . . , or those labelling the additional 
functions. 

From now on we will treat the additional functions Ja (which may or may not be given 
together with some explicit functional dependence among themselves and the /^j/a) as ad- 
ditional fundamental differential invariants, and unless explicitly mentioned, "fundamental 
invariants" means both I^^x and Ja- 

Remark 1. It is perfectly possible to construct from coordinates. But once these 
forms have been constructed, we should treat them as fundamental and forget the coordin- 
ates. All geometrical properties of the system will be derived through the invariants I^uX- 
This is the whole point of the method of moving frame. ^ 

Remark 2. We can take the indices i,j,k... to run over all values of //, z^, A . . . , and 
then the assumption of the existence of this subset of indices will be valid for all systems. 
The utility of our assumption, on the other hand, will be clear in a moment. ^ 
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Remark 3. In cases where we have additional functions Jq,, these are considered further 
constraints: two systems are considered equivalent if and only if we can set uj^ = oj'^ and 
Ja = J a- The incorporation of the possibility of additional function into our algorithm 
gives us great flexibility in application of our theory. For example, if for example, we have 
a system with one- forms oj^ spanning the cotangent space but are not independent, we can 
choose a subset of these one forms as the co-frame and the rest can be written in terms 
of the CO- frame with the help of additional functions, and thus this system can be treated 
by our algorithm. On the other hand, if a system has certain one-forms uj^ which do 
not span the cotangent space, we can choose arbitrarily some other one-forms in order to 
form a coframe, and then prolong the problem to obtain a coframe on a suitable principal 
bundle, [23] has more details on this. Thus in this case our algorithm is also applicable. ^ 

2.2 The derived invariants 

We can differentiate the fundamental invariants: 

where the functions If_iux-s are called first order derived invariants. As uj^ form a co-frame, 
they are uniquely determined in terms of the co-frame. Derivations can be carried out 
further: 

where the functions I^,yx-Sy are called second order derived invariants. Derived invariants 
of all order are defined recursively. 

By our assumption, the one-forms lv^ other than oji are labelled using indices i,j,k . . . , 
and we denote such forms hy uji,ujj, . . . where I,J . . . are multi-indices in i,j,k ... . Thus, 
for any fundamental or derived invariant Ii-j, we have 

dll;j = Il;Jk^k + Ir,J;Ki^K, 

which is nothing more than the above formulae rewritten using different indices. We assume 
that the derived invariants Ii-j-k can be expressed linearly in terms of Ii-jk and lower order 
invariants. 

Remark 4. As before, if we take i,j,k... to run over all possible indices, then the above 
assumption is trivially verified since there is then no Ii-j-k- In almost all applications, 
however, the uji will be the horizontal one-forms in a principal bundle and the above 
assumption means that only the horizontal derivatives can be independent, i.e., we have a 
connection at work here. The derived quantities Ii-jk are then the covariant derivatives of 
I J- J in the bundle if the bundle is reductive [26]. ^ 

2.3 Algebraic relations 

Next we need to take into account algebraic relations of the invariants. When we write down 
the structural equations, the fundamental invariants appearing in the structural equations 
are subject to certain relations (symmetries). For example, for the most general structural 
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equation (|2.1|) . the symmetry is /^^.a = —IiiXv When we stipulate the additional func- 
tions, they may also be subject to certain relations. We call such algebraic relations the 
fundamental algebraic relations. 

We can derive the structural equations: 

= (fuj^ = dlf,u\ Auji, Au;x + I^vxduj^ A dujx - I^j^yxu^ A dujx, 

and, after using the structural equations themselves and the resolution of dl^^x in terms 
of derived invariants, we obtain 

= F^yx^^ Aujy AuJx, 
where F^^x are functions in the differential invariants. We call the algebraic relations 

the Bianchi algebraic relations. 

If Ij-j is a differential invariant, deriving it we get 

dll;j = Il;Jk'^k + Ir,J;K^^K, 

and deriving again we get 

= Il-Jkl^^k A W/ + Cijkl^k A W/ + AijAk^^A A Wfe + BijAB^^A A LOb, 

where the functions Cjjki contains only invariants of lower order. We make the further 
assumption that Ajj^k = and Bijab = identically. 

Remark 5. Again, if i,j,k... runs over all values the assumption is trivially verified. 
If, as we have mentioned, the choice of indices i,j,k... comes from the existence of a 
connection, this assumption can also be verified easily. % 

The relations 

Ir,jki = hjik + ClJkl 

are called the generic algebraic relations. 

Obviously, if an invariant // is actually one of the additional functions in the system, 
there will not be any Bianchi algebraic relation for it. 

If i? = is an algebraic relation, we can derive it to obtain 

dR = RiUJi + Ra^a = 0, 

and we have the new relations 

i?. = 0, Ra = Q 
which are called the derived algebraic relations. 

Remark 6. As before, the relation Ra = is usually either vacuous or an identity. ^ 

The defining relations, the Bianchi relations, the generic relations and their derived 
relations are all the algebraic relations of the invariants. 
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2.4 The involutive seeds and the degree of arbitrariness 

First we need to define the concept of covering of invariants. In the following, when we say 
that a function / is expressible using the functions ga, a € I where / is some indexing set, 
we mean that we can find an explicit function / such that 

/ = fi9ai,ga2,---,9ak)^ ttl, Q!2, . . . , ajfc G /. 

In particular, if / is a constant, then it is expressible using the empty set. 

Let S" be a set of invariants (the covering set) and T be another set of invariants (the 
set to be covered). From 5, we form the set S' by 

1° adding all elements expressible using the elements of S to S' 

2° adding all elements Ij-j such that Ii-jk for k can be expressed in terms of elements 
in S, 

and we repeat this procedure, using S' in place of S, to obtain S", etc. This procedure 
may eventually stabilise. In practice, due to the way we use the set S^''\ we only need to 
add terms that are independent: sec below. 

Suppose at a certain stage we have the set S^'^K We remove all elements of T that can 
be expressed in terms of elements in S^'^^ to form a set T^*^). If every element Jj-k of T^''^ 
are such that Jj-kl for all multi- indices L of a certain length / are expressible in terms 
of the empty set, i.e., are constants, then S covers T. In particular, S covers T if T^*^) is 
empty. 

This definition may seem a little long-winded and peculiar, especially the part dealing 
with non-empty T^^). As a result of this definition, if a set of invariants Ii-j has Ii-jl = 
constant for all multi-indices L of a certain length, then it is covered by the empty set. 

A system of involutive seeds is a set of differential invariants that satisfies the following 
conditions: 

11. (Covering.) All invariants occurring in the structural equations = • • • and the 
derived structural equations do;^ = • • • , as well as in any algebraic relations to be 
enforced, are covered; 

12. (Independence.) All invariants in the set are algebraically independent, namely 
there is no non-trivial relations among them, possibly also involving invariants 
that are covered by the set of involutive seeds, that is derivable from the algebraic 
relations that we mentioned. Hence, no proper subset of this set is a covering of 
the above-mentioned invariants. 

13. (Derivation index.) Every invariant in the set contains at least one derivation index. 
(This condition can be non-vacuous only for invariants coming from additional 
functions.) 

Now we will give the indices i,j,k . . . ,n an ordering labelled by the first few natural 
numbers, i.e., a function s such that s{i) G N, which amounts to a permutation, or re- 
labelling, of the numerical indices if the indices themselves run from 1 to some number 
without gaps. An ordering will be called an involutive ordering if the following two condi- 
tions are satisfied: 

01. (Maximality of lower indices.) If two invariants In and Jjj occur in a relation for 
some value of i and j and i > j, then In must not be an involutive seed. 
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02. (Counting condition.) For every involutive seed Iijk...i, for all indices m < I using 
this ordering, the invariants Iijk...m are also involutive seeds. (Note: they must 
occur explicitly, not merely expressible as linear combinations of others.) 

The names given to the conditions appearing in parentheses above will become clear when 
we prove our algorithm. 

Remark 7. That an ordering we write down is an involutive ordering with respect to 
a system of involutive seeds is the test of our algorithm, and we will obtain information 
about the degree of arbitrariness only if this test passes. In the next section we will show, 
by examples, of how to proceed in order to have a good chance of arriving at an involutive 
ordering. ^ 

We are now nearly done: armed with an involutive ordering for a system of involutive 
seeds, in many cases we can read off the degree of arbitrariness of the system directly. Let 
us describe one further test: 

R. (Rank condition.) We can find a number of algebraically independent invariants Iq, 
a = 1, 2, ... , which are not additional functions introduced and which when derived 
give 

and the rank of the matrix C^^ is equal to the number of one- forms oj^. 

Assume that this test is satisfied. Let S' be the subset of S such that all elements are such 
that their last index is maximal with respect to the involutive ordering among elements 
of the sets. Let this index be d and let k be the number of elements of S'. Then we say 
that the system has degree of arbitrariness k, occurring at dimension s{d). If the set of 
involutive seeds is the empty set (this is possible, since we have seen that a set can be 
covered by the empty set), then the system does not have any degree of arbitrariness. It 
may or may not be inconsistent. 

Remark 8. Another way of phrasing the same result is that the system depends on k 
independent, arbitrary functions of s{d) variables. See section HI ^ 

Granted, there are systems for which the test R fails. For such systems, we can still 
apply our algorithms, but the number we obtained is only meaningful if the system does not 
involve additional functions that we put in by hand, and still the number is only guaranteed 
to be an upper bound of the degree of arbitrariness of the system. 

Remark 9. As we will see, if the test R fails, then the question we are posing are not 
formulated on the best space possible: we can reduce the problem into a lower dimensional 
one, for which the test R holds and we can obtain the degree of arbitrariness, not merely 
an upper bound. If R fails and we have additional functions defined on the space, we may 
be able to use the additional functions to specialise the moving frames so as to reduce 
the dimension. Note that such a reduction of dimension need not be a reduction of the 
dimension of the base manifold: for example, in the principal bundle over a manifold 
Riemannian, we may be able to reduce the bundle from an SO{n) bundle to a suitable G 
bundle where G is a subgroup of SO{n), and in extreme cases G may even be a discrete 
group. 
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As any well-formulated problem will satisfy the test R, most often the verification of 
this condition is easy. ^ 



2.5 Adding "equations of motion" 

After we have obtained the degree of arbitrariness coming from the geometry of a system, we 
have effectively solved parts of the so-called Cauchy problem for the system: the involutive 
seeds that have their last index maximal are those that wc can specify arbitrarily, which, 
together with suitable boundary conditions, will determine the system uniquely. 

However, these functions that we can arbitrarily specify are not the only possible choice, 
and often they are not the functions that we really want to specify: most of the time they 
are derivatives of too high an order. Usually what we really want to do is to specify some 
functions of the low order invariants. The question is then: can we specify these functions 
arbitrarily? 

In a guise, this is precisely the problem of adding equations of motion to a geometrical 

(kinematical) system. The equations of motion are usually differential equations, but as 
higher order invariants are nothing more than derivatives of the fundamental equations, 
these equations of motion are, in our framework, just additional algebraic relations. Hence, 
by taking these equations as additional defining relations for the differential invariants, 
our method is still applicable: the new degree of arbitrariness, after the introduction of 
equations of motion, signify the Cauchy data for the equations of motion. 

Note that after the introduction of additional relations the old system of involutive 
seeds and ordering may not apply at all: for example, a completely different ordering may 
be needed. The new system may also fail to admit any system of involuliA-c seeds and 
ordering: this is in particTilar the case if the eqTiations of motion introduced is potentially 
contradictory. These will all be illustrated by examples later. 

Remark 10. Real, physical equations of motion have additional nice properties: they 
always reduce the dimension at which the degree of arbitrariness occurs by exactly one. 
See the examples later. % 

2.6 The directions of evolution 

There is one additional bit of information that we can obtain from our manipulation of 
indices. Assume that for a system, our algorithm applies with the condition R satisfied. 
Then for a system of involutive seeds together with its involutive ordering, let 5* be the set 
of last indices of the involutive seeds, and let T be the set of all possible indices. Then the 
set of indices T — S gives the directions of evolution of the problem. The interpretation is 
as follows: if we consider the Cauchy problem of the system, then the system is specified 
completely by specifying d functions on a A; dimensional submanifold of the manifold in 
which the system is defined, where d is the degree of arbitrariness and k is the dimension 
at which the degree of arbitrariness occurs. But in general, this submanifold cannot be 
chosen at will: it must be transverse to the system of vectors having indices taken from 
the set T — S for a certain choice of involutive seeds and ordering (this choice is in general 
not unique). As a consequence of our assumptions that the one-forms splits into two sets 
LOi and LJa and all independent invariants take the indices of Wj, we see that all directions 
corresponding to Ua are automatically directions of evolution. If the space for which the 
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system is defined is a principal bundle and uji are the horizontal forms whereas uo^ are the 
vertical forms, this just affirms the fact that the data in the bundle is completely determined 
once we specify the data on a section of the bundle. 

There is a further constraint on the choice of the submanifold: it must not contain the 
so-called characteristic directions. For this constraint, see the proof of our algorithm. 

Remark 11. For real, physical equations of motion, the single additional direction of 
evolution when compared with the system without the equations added must be the time 
direction. ^ 



2.7 A little adjustment in case of failure 

There are cases for which our algorithm is not applicable, since we will not be able to find 
any system of involutive seeds and involutive orderings. Most often, the problem is of the 
following kind: wc have two invariants Ii-^jk and Ij/.j//^., and for k = a, Ii-ja is independent 
and I'li.jig^ is not, whereas for k = b, is independent whereas Ij-jb is not. Thus we 

can take neither a > b nor a < b for our involutive ordering. 

In most such case, the solution is easy: we just take the forms cja and cjf, and replace 
them with u'^ = acoa + P^b and = ^yua + Sub, with 

and it is easy to show that, for this new co-frame, which is otherwise completely equivalent 
to our old one, we can either take Ii-^ja and Ij/.j/^ to be independent, or Ij-ji, and Ip.jif,- 
Of course, for complicated systems involving a large number of offending one-forms, this 
might get quite complicated. 



3 Examples 

If the description of the algorithm above seems too abstract, the examples below will show 
how easy it is to apply it and how little the work required is. 

3.1 Riemannian geometry and Einstein gravity 

A Riemannian geometry is defined locally by specifying a non-degenerate, positive definite 
symmetric bilinear form. (The assumption of positive definiteness is not required for what 
follows and is added only to simplify the exposition. However, there will be differences 
when we discuss the equations of motion: see the remarks later.) In dimension n, such a 
bilinear form has n(n + l)/2 independent components. By a not very precise argument, 
since two Riemannian metrics g and g' on M and M' are considered equivalent if there is 
a function f : M ^ M' such that f*{g') = g, the real degree of arbitrariness is the number 
of independent components less the number of coordinates, namely n{n — l)/2. 

In the moving frame formulation, a Riemannian geometry on M is locally defined by a 
principal bundle tt : P — M, which is locally the same as SO{n) x M. A suitable coframe 
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on P is given by the system of horizontal forms uoi and the vertical forms ujij, with oj^j 
forming a representation of the Lie algebra so(n). The structural equations are 



(3.1) 



duji = —ojij A Uj, 

dcjij = —ujik A Ukj + ^RijkiuJk A uii, 



which can also be written in the matrix form: 

dio = —Lo A w + ri, 

where 

oj=(^ 0\ ^^/O 

\0Ji UijJ ' \0 RijkiUJk A uji 

making the Lie algebra structure explicit. The fundamental invariants are Rijki, with 
defining symmetries 

(3.2) R'ijkl = —Rjikli Rijkl = —Rijlk, 

the first coming from the fact that coij is a representation of so(n), and the second due to 
how Rijkl appears in the structural relation. Obviously we want to take Ui as the basic 
forms, as everything else is labelled using their indices. The definition of the first derived 
invariants 

(3.3) dRijkl ~ Rijkl-jmJ-^m Rmjkl^im Rimkl^jrn Rijuil^kni ~ Rijkm^lm 

shows that our first assumption about the basic forms is satisfied, and Rijki-m are really the 
covariant derivatives of Rijki- The definition of higher order invariants are straightforward. 
For the Bianchi relations, d^oji = gives 

(3.4) Riijki] = 0, 
whereas d^ojij = gives 

(3-5) Rij[kl;m] = Oj 

i.e., they are just the usual Bianchi identities. The generic relations and the derived rela- 
tions are obvious. 

For this system, the invariants Rijki-m cover all the invariants occurring in the struc- 
tural equations and the derived structural equations. Our task now is to find a system of 
involutive seeds and an involutive ordering. For this system, no index is special and hence 
we take any ordering in which i,j,k... run from 1 to n. 

It is best to arrange the set of invariants themselves systematically in order to arrive 
at a system of involutive seeds. Since all derived algebraic relations must be taken account 
of, it is best to arrange the invariants order by order. We will use the following strategy: 
when using an algebraic relation relating invariants, we prefer to take the invariants having 
the following property as the independent ones: 

1° They are of the lowest order; 

2° Their indices are of decreasing order, if possible. 
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This "if possible" is intentionally vague since it is not an requirement, only a guideline 
aiming to make the condition 01 hold. For the defining relations ()3.2p . this suggests that 
we take the following to be independent with respect to this relation: 

Rijkl such that i > j, k > I. 

For the first Bianchi relation (|3.4p , we will use a now well-known trick to make our work 
simpler. We expand the identity 

Ri[jkl] — Rj[ikl] - Rk[ijl] + Rl[ijk] = 

and use ()3.2p : this gives 



Rijkl Rklij ) 

SO in addition, we can require that the independent terms satisfy 

Rijkl such that i > k. 

However, this does not use up all the relations Ri[jki] = 0. Note that unless all four indices 
are different, this relation is trivial. Assuming all indices different, and i > j > k > I, 

Rijkl — Rikjl + Riljk = 

contains two terms for which the second index is greater than the fourth and one term 
for which the second index is less than the fourth. Hence we can also require that the 
independent terms satisfy 

Rijkl such that j > I- 

There are now no more relations on Rijki- Incidentally, this procedure counts the number 
of algebraically independent terms of RijkV- 

(3.6) i> 3, k>l, i>k, j>l 

giving 



4J \3J \2J 12 
independent terms. 

Next we come to the first order invariants, Rijki;m- First we need to take care of 
the derived relations. This means that we may take Rijki-m to be independent only if it 
results from the derivative of an independent fundamental invariant, in other words, we 
still require (j3.6p . The only additional relation that these invariants are subject to are the 
second Bianchi relations ()3.5p . Using a completely analogous reasoning as the one we used 
for the Bianchi identity, we require that independent terms satisfies in addition /c > m, in 
other words, the requirement on the indices is now 

(3.7) ^ > J) k > I, i > k, j > l, k > m. 

There are no generic algebraic relations since we only have one derivation index. 

Having obtained all the algebraic relations we need, we can now verify the condition 
R: the expansion (j3.3p shows that we really have too many invariants for it to fail. For 



14 



example, for the invariants Rn^j]^, Rninli Rn,n-i,n,n-2-, the matrix in question is already of 
maximal rank: it suffices to pay attention to the terms linear in Uij in the expansion of 
d^ikjk' terms linear in tOi in the expansion of dR^^f^i, and the term linear in cOn in the 

expansion of dRn,n-l,n,n-2- 

Now we take the independent terms of Rijki-rn the involutivc seeds: they obviously 
cover all of Rijki ;m since they are obtained by applying all the algebraic relations on Rijki;rm 
and since all of Rijki;m are covered, all of Rijki are as well. The involutive ordering is the 
ordering that we have been using all along, and that this is an involutivc ordering is obvious 
from the fact that the only requirement on the last index m is A; > m, so if ^? < m, we also 
have k > p. 

Now, what is the maximal value that can be taken by m? Let us try n, which is the 
dimension of the manifold. Since we require 

i > k > m = n, 

we have 

i = k = m = n, n > j > I, 

so the number of involutive seeds with m = n is the same as that of a symmetric bilinear 
form of dimension n— 1, namely n(n— 1)/2. This is the answer that we have been expecting. 

Note that we can take the independent invariants at any order to be the involutive 
seeds. If we do this, besides the derived relations, we need to uphold the generic relations, 
which means that for our purpose the derivation indices must be non-increasing. Hence, 
for example, at the third order, the seeds with the maximal last index are Rninj;nnn, with 

n> i>j, 

giving the same answer n(n — l)/2. 

Remark 12. One is justified to have the feeling that for this problem, it is not necessary 
to formulate the problem in the language of moving frames at all. Indeed, let Rfj,i,p\ be the 
Riemann tensor in the usual language of tensor analysis, then our fundamental invariants 
are just 

Rijkl — f^i ^j^k^l R^ivpXi 

where is the linear transformation that transforms the coordinate frame into an or- 
thonormal frame (they are the "vierbein"). And for our purpose, the symmetries of Riju 
and those of R^upX are equivalent since we ignore all non-indcpcndcnt terms. 

It is indeed true that, provided one is very careful, the algorithms can be applied to 
tensors. But using tensor analysis, being careful is extremely difficult: with the moving 
frame, it is simple to exhaust all the algebraic relations applicable to the system, and 
on the other hand, using tensor analysis it is very easy to overlook one or two relations. 
Structurally speaking, the methods of moving frame are much more rigid than the methods 
of tensor analysis, and there are crazy things that can be done to tensors using tensor 
analysis but are forbidden using moving frames. Hence, to ensure the applicability of the 
algorithm to general tensor systems, many more tests and conditions are required, which 
is far more difficult than rewriting the system in the language of moving frames. ^ 
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Riemannian geometry with torsion. If, instead of writing the structural relations (|3.1|) , 
we write 



(3.8) 



duji = —LOij A LOj + ^TijkOJj A OJk, 

dcOij = —OJik A UJkj + \Rijkl^k A CJ;, 



we have a Riemannian geometry with torsion T^-fc = —Ti^j- We can carry out the above 
procedure as before: the algebraic relations for Rijki are, up to terms in Tijk and Tijk-^i, 
which has no effect in our theory, exactly the same as before, and there are no separate 
Bianchi relations for Tijk- Thus, for the involutive seeds, we use Rijki-m with symmetry as 
before, and Tij^-i. As Tij^-i satisfies no Bianchi relation of its own, there is no restriction 
on and for I = n the number of seeds is just the number of independent Tijk. As j and 
k are antisymmetric by the defining relation, the additional contribution is n^(n — l)/2. 

Actually, since we already know the degree of arbitrariness of Riemannian geometry 
without torsion, we can reason as follows: the structural equation ()3.8p defines the two 
one- forms and uJij, from which we can define uniquely the two one- forms LOi and uj'-, 
such that these two forms satisfy (|3.ip . by setting 

Hence the system (|3.8p is completely equivalent to the system p.ip together with the 
additional function Tijk with the relations 

Tijk — Tikj, 

from which we can easily deduce the additional contribution to the degree of arbitrariness. 



ij + { — Tikj + Tjki + Tkji)Ldk- 



Ricci-flat spacetime. Einstein equations. Now let us return to the theory without 
torsion, and focus on Riemannian spaces that are Ricci-flat. Such spaces are defined by 
the constraint 

k 

Note that this constraint is equivalent to the vacuum Einstein equations 

Rij ~ \^ijR = R = 

i 

and hence Ricci-flat spacetimes are just those satisfying Einstein equations in vacuum (and 
without cosmological constant). (We have pretended that Einstein's equations are defined 
using a positive-definite metric. The construction and analysis, done properly using metric 
with the correct signature, is essentially identical to what we do here.) Remark that these 
constraints, when written in terms of the metric, are differential equations, but for us they 
involve no differentiation at all. 

For this system, we can make the following changes to our calculation for the general 
Riemannian case to obtain the new degree of arbitrariness: for i,j < n, the equations we 
have added are 

Rlilj + R2i2j + • • • + Rninj = 0, 
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so now, by the condition 01, Rninj for i,j < n are no longer considered independent terms. 

For i = n, j < n — 1, we have 

Rlnlj + R2n2j H 1" Rn-l,n,n-l,j = 0, 

SO Rn,n-i,n-i,j for j < n — 1 are no longer considered independent. For i = n, j = n — 1, 

Rl,n,l,n-1 + R2,n,2,n-1 H + Rn-2,n,n-2,n-l = 0, 

so Rn,n-2,n-i,n-2 is no longer considered independent. Finally, for i = n, j = n, 

Rnlnl + Rn2n2 + • • • + Rn,n-l,n,n-l = 0- 

However, all of these terms have already been declared dependent in the case of i,j < n. 
Hence we need to substitute them with independent terms. After this substitution, we see 
that Rn-i^n-2,n-i,n-2 is no longer independent. 

In summary, the independent (normal) components of the Riemann tensor are now Rijki 
for which 

i > j, k > I, i > k, j > I, 
from which we exclude the following 

Rninj 1 Rn,n—l,n— 1,1-1 Rn,n—2,n—l,n—2j Rn—l,n—2,n—l,n—2- 

We also need to take the equations dRij = into account. The derived equations 
only tell us that independent derived invariants comes from derivation of independent 
fundamental invariants. 

Now let us apply our algorithms. As in the unconstrained case, the verification of 
condition R is easy. The largest possible last index for the seeds is now n — 1, and we have 
the following possibilities: 

Rn—l,i,n—l,j;n—l} Rn,i,n—l,j;n—l- 

The first possibility gives 

(n-l)(n-2) ^ 
2 

terms (we need to exclude Rn-i^n-2,n-i,n-2-n-i>, and the second possibility also gives 

(n-l)(n-2) ^ 
2 

terms (we must have i < n — 1, and exclude Rn,n-2,n-i,n-2)- Hence the degree of arbitrar- 
iness is now 

n{n — 3) 

at dimension n — 1. 

For the general Einstein equation 

Rij 2 ^'^j ^ ^^3 ~ J ' 
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the degree of arbitrariness in the general case is the same as in the vacuum case, as can be 
verified by a similar but slightly more complicated calculation, by using the expression for 
the Ricci scalar R. 

It should be noted that n{n — 3) is always an even number, and this shows that the 
Einstein equations are well-behaving dynamical equations suitable for describing particles: 
the number n(n — 3) /2 is of course just the degree of freedom of the graviton, which when 
n — 4 yields the familiar answer 2. 

Let us now briefly consider Einstein-Cartan theory, in which we have Riemannian geo- 
metry with torsion. In addition to the usual Einstein equations, we couple the torsion to 
spin density directly: 

so that all degree of arbitrariness generated by Tjjfc are killed. Thus, the theory has exactly 
the same degree of arbitrariness as the usual Einstein theory. Physically, we say that "spin 
density does not propagate" : there are no "spin waves" . 

Remark 13. As treated here, the direction of evolution of Einstein's equations can be 
taken to be any direction, i.e., any submanifold can be used as a Cauchy surface. But this is 
only due to the fact that we are formulating the theory as if the metric is positive-definite: 
in this case the system admits no characteristic directions. If we use the metric of the 
correct signature, or what amounts to the same thing, if we write the structural equations 
(j3.1|) as the structural equations for the group SO{n — 1,1) x M instead of SO{n) x M, 
which in practice just adds a few minus signs here and there, one important change to 
the theory is that the theory now admits characteristic directions, which are just the null 
directions. As we have discussed, the Cauchy surface must not contain characteristics, and 
hence the signature forces the time direction to be the direction of evolution. ^ 

Remark 14. If we compare the theory before and after adding the Einstein equations as 
constraints, we see that the dimension at which the degree of arbitrariness occurs is exactly 
one lower than the unconstrained theory. This property, together with the property men- 
tioned above about the time direction, are the requirements for any physically meaningful 
equations of motion: the theory must be neither "over-determined" nor "under-determined" 
by a suitable Cauchy data. ^ 

Historical Note. For n > 4, the degree of arbitrariness n(n— 3) we obtained for Einstein's 
equations is greater than the number of degree of arbitrariness for a n — 1 dimensional 
Riemannian space, but less than the degree of freedom of a n — 1 dimensional Riemannian 
space plus a symmetric bilinear form (For n = 3, the degree of arbitrariness by our formula 
is zero: we can see this even without calculating since Rij = implies Rijki = when 
n = 3). Hence we may speculate that the degree of arbitrariness that we have obtained 
decomposes into the metric on the n — 1 dimensional "section" and the second fundamental 
form on this section, subject to certain constraints. This suggests that we can try to use 
a metric on a hypersurface and a second fundamental form subject to certain constraints 
as the Cauchy data for the Einstein equations. That such a specification actually works is 
shown by Choquet-Bruhat [13] in her study of the Cauchy problem of Einstein's equations: 
for a modern account, see ^14j. However, since the second fundamental form is subject to the 
Gauss-Codazzi equation, and in Choquet-Bruhat 's approach there are further constraints 
on them (the momentum constraint and the Hamiltonian constraint), it is impossible to 
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derive the degree of arbitrariness in this way without a large amount of further work. 
On the other hand, before the work of Choquet-Bruhat, Cartan has aheady derived the 
correct degree of freedom of the vacuum Einstein equation n(n — 3) at dimension n — 1 by 
a rather ad-hoc reasoning, which is comphcated by the arbitrariness of coordinates, in his 
correspondence with Einstein: see [15], letter xxii. ^ 

3.2 Gauge theories and Yang-Mills equations 

Here we want to find the degree of arbitrariness of a classical gauge theory over a Rieman- 
nian geometry. A gauge theory on a Riemannian manifold is usually specified by writing 
down some gauge potentials Aai, for which a is a group index (omitted if the group is one- 
dimensional) and z is a spacetime index, and Aai depends on the spacetime coordinates 
only. This requires pn functions to define, where p is the dimension of the group and n is 
the dimension of spacetime. However, there is also the "gauge invariance" of Aai that must 
be taken into account: thus, the true degree of arbitrariness is, by this intuitive argument, 
p{n — 1). 

Let us calculate this number by our algorithm. In the formulation of moving frames, a 
gauge theory is constructed as follows: we have the structural equations for the Riemannian 
space (|3.ip . and we couple a Lie group G to it. Let us assume that the Lie group has 
Maurer-Cartan structural equations 

dOa = —CabcOib A Oc, 

where aa are the Maurer-Cartan forms and Cabc are the structure constants for the Lie 
group. We form the product space of the Riemannian principal bundle and the Lie group, 
and change the structural equation to 

d'ya = -Cabclb A 7c + ^FaijUJi A LOj. 

The defining relations for Faij are as follows: for the index a, the symmetry is the same as 
that of the form Ua- For the index z, j, we have Faij = —Faji- There is also an additional 
relation occurring at first order, namely 

dFaij = Faij-k^k + Gujlb + (terms in tOij), 

where G^j are functions of Faij-, the exact form depending on the group. For example, for 
the group SO{p), the additional one-forms are jab with a, 6 antisymmetric, 

dFabij = Fabij;k^k F^cbijlac Pacijibc F^abkj^ik ^abik^jk- 

Basically, the form of these relations just means that the group directions are not dynamical, 
as differentiations in these directions do not generate new independent invariants. 
The Bianchi relations for Faij is 

■'^ a{ij\k\ ^1 

and the rest of the relations are obvious. Thus, we can take the independent terms of Faij 
to be those that have i > j, and using the Bianchi relation, take the independent terms of 
Faij;k to be thosc that have 

i > j, i > k. 
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Then we take the set of involutive seeds be the independent Rijkim ^ 

before, and also 

the independent Faij-k- The involutive order is the order that all indices a,b,c... are 
considered greater than i,j,k . . . . The number of terms Fani-n is thus 

p{n - 1), 

where p is the dimension of the Lie group, since the index a in Faij is a group index. The 
condition R can be easily verified as before. Hence this is the additional contribution to 
the degree of arbitrariness, occurring at dimension n, as the not-so-precise argument at the 
beginning of this section shows. For example, for SO{2) gauge theory (electromagnetism) 
in dimension 3 -|- 1, this number is simply 3. 

Yang— Mills equations. We now study the degree of arbitrariness of adding to the above 
system the classical Yang-Mills equation, which in our notation, reads 

Faij-i = source terms, 

i 

which has a total of pn{n — l)/2 equations: more than the degree of arbitrariness of the 
original system. For the equations of motion, when j < n, the constraints are just that 

F 

are no longer independent by condition 01. When j = n, we have (omitting the group 
indices) 

Fln;l + -^2n;2 + • • • + Fn-l,n;n-l = ' ' ' , 

SO Fn^n-i;n-i IS no longer independent. Under Yang-Mills equation for the gauge fields and 
Einstein's equations for the gravitational field, the additional contribution to the degree of 
arbitrariness from the gauge fields is therefore the number of remaining normal terms of 
Faij;n~i, which are 

Fa,n,j;n^l (j = 1 , • • • , n - 2) , Fa,n-l,j;n-l (j = 1, . . . , n - 2), 

giving the degree of arbitrariness 

2(n - 2)r. 

occurring at dimension n — 1. This gives the number of fields we must specify on a hypersur- 
face to have a well-defined Cauchy problem for Yang-Mills equations coupled to Einstein 
gravity, or in flat spacetime. Note that it is essential to impose the Einstein equation as 
well, otherwise the 2(n — 2)r degree of arbitrariness we get here at dimension n — 1 will be 
eclipsed by the degree of arbitrariness of the Riemannian geometry at dimension n: more 
on this in the next section. Again, an essential feature is that when going up one dimen- 
sion, we require two additional copies of the Lie algebra: this corresponds to one degree of 
freedom for the boson. When r = 1 and n = 4, we see that the photon has 2 degrees of 
freedom. 

Historical Note. For specific dimensions and specific groups, Estabrooks [T7] has set 
up explicit exterior differential systems for the Yang-Mills equations using coordinates, 
and, with the help of computer algebra programs, calculates the Cartan characters for the 
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systems (which, as we shall see in the proof of our algorithm, is intimately related to the 
degree of freedom). For example, for SU{2) Yang-Mills equations of dimensions 3,4,5,6, 
the last non-vanishing Cartan characters occur at dimensions 2, 3, 4, 5, and are 9, 15, 21, 27. 
However, since Estabrooks used coordinates, the gauge degree of arbitrariness (3 for SU(2) 
Yang-Mills theory) is still present. If we subtract it from his answers, we get 6, 12, 18, 24, 
which is just our answer. The Cartan characters obtained for Maxwell theory for dimensions 
3, 4, 5, 6 by Estabrooks can also be shown to be in complete agreement with our result by 
analogous reasoning. Of course, our algorithm is so simple so that it is unnecessary to 
resort to computers for the calculations, we not need to reason with the gauge degrees of 
arbitrariness since no coordinate is used, and our result does not depend on the dimension 
nor on the geometry of the underlying space (Estabrooks considered only theories set up 
in fiat spaces). % 

3.3 Scalar field theory on Riemannian manifold 

We have checked that our algorithm gives the correct answers for gravitons and gauge 
bosons. Let us now very briefly check the case of scalar fields. In our approach, adding 
a scalar field corresponds to having a free function to a Riemannian manifold. From our 
discussion about additional functions, it is obvious that at the kinematical level, the extra 
degree of arbitrariness is 1: denoting the field be the scalar function /, this degree of 
arbitrariness comes from f-n- Now we study what happens when we specify the equations 
of motion for the scalar field, i.e., its dynamics. 

We add the Klein-Gordon equation, which is, in moving frames, the single equation 

n 
i=l 

Assuming that we have already killed all degrees of arbitrariness coming from the Rieman- 
nian metric. For any system of involutive seeds and ordering, f-nn is no longer considered 
a seed. Hence, by the generic relations, the degree of arbitrariness on a Cauchy surface 
comes from the two terms 

/;n,n— 1) f;n—l,n—l 

and thus the degree of arbitrariness of a scalar field under Klein-Gordon equation is exactly 
2, independent of the dimension of the manifold, occurring at dimension n — 1. Of course, 
this translates to a degree of freedom of scalar particles | = 1- 

3.4 Newtonian rigid motion 

Of course, the point of this algorithm is not for studying the degrees of freedom of element- 
ary particles, which arc already known. Thus let us now consider constrained systems, the 
first example of which being the familiar Newtonian rigid motion. 

Intuition tells us that a Newtonian rigid motion has six degree of arbitrariness. As a 
Newtonian rigid motion is locally just a rigid flow (i.e., shear-free and expansion-free flow), 
this is also the degree of arbitrariness of such a flow in fluid mechanics, but note that in our 
language, this degree of arbitrariness occurs at dimension 1, even in the fluid description. 
Below we shall derive this rigorously using our algorithm. This example also shows a few 
pitfalls when applying our procedure. 
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Newtonian motion is formulated on a space with a Galilean connection defined on it. 
The connection matrix can be written 

/O \ 
T 

\9i U!i UJijJ 

with the additional constraints 

dT = 

which guarantees the existence of absolute time, 

diOij = —Ldik A Ld^j 

which guarantees the existence of flat space, and 

{doJi + bJij A UJj) A 9i = 

which guarantees the absence of velocity-dependent gravitational effect. For how these 
equations are derived, see [IT], or the English translation |12j. In the usual coordinates, 
this matrix is of the form 

0\ 

dt , 
dxi Gidt 0/ 

but we will not use any of such coordinates. This moving frame can be adapted to the 
rigid motion of our interest if we align the co-frame r with the direction of motion. In the 
language of fluid mechanics, the flowlines, characterised by r, drags the spatial structure 
along while keeping them unchanged. The structural equation of a rigid motion is 

' dT = 0, 

dOi = —uji At — uji-j A 6-j, 
< 1 

doJi = —ojij A UJj + TijOj At + ■^Tiji^.Oj A 6^, 

, dujij = —LOik AuJkj, 

which are just the structural equations for Newtonian spacetime (the one- form r implicitly 
has the index 0), together with the decomposition 

LOi = KiT + MijOj. 

For more details, see [21]. Thus, properly speaking, the co-frame is really formed by r, 6j 
and LOij'. the Ui are really just convenient one-forms to work with. Roughly speaking, Ki 
signifies the acceleration and Mij the rotation. The fundamental invariants are hence Fjj 
and Fjjfc, with a defining relation 

together with Ki and Mjj, with Mij = —Mji. The condition of the absence of velocity- 
dependent gravitational effect implies 

^[ijk] = 0, = 0. 
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Here some of the Bianchi identities are obtained by expanding 

duji = d{KiT + Mij0j) 
using the structural equations, which gives 

(3.9) Mij-k = 0, Ki,j = Mij,o + MikMkj. 

The other Bianchi relations are easily derived to be 

^ijk = 0, ^i[j;k] = 0, 

but as = 0, rjj.fc is totally symmetric in all three indices. It is easy to check that we 
again have an abundance of invariants, so that condition R is satisfied. 
It is attempting to take the set 

Ki-j, Ki-o, Tij-o {i > j) Tij-k {i>j> k) 

as the set of involutive seeds. However, this does not work: the second relation in ()3.9p 
means that the condition 01 fails for Ki-j if the involutive ordering has as the smallest 
index, and the condition 02 fails for Tij.^ and Tij-Q if is not taken as the smallest index. 
Instead, we must take the system of involutive seeds to be 

Mij-o, Ki-o, Tij-o {i > j) ^ij;k {i>j> k) 

and specify the involutive ordering as < i,j,k,..., the ordering among i,j,k... being 
immaterial. We see that the degree of arbitrariness is 1, given by the term Ti^—\ Yi—i;n—i^ 
occurring at dimension n. 

But actually this system is not exactly what we really want to use if we want to calculate 
the degree of arbitrariness of Newtonian rigid motion: the presence of the derivatives of 
Tij means that the degree of arbitrariness of the gravitational field is still present. The 
degree of arbitrariness of the gravitational field is also 1 at dimension n: we can calculate 
this separately but it is most easily seen if in the above system of involutive seeds we 
simply remove Mij-o ^ind Ki-o. We say that the degree of arbitrariness of the rigid motion 
is completely eclipsed by the degree of arbitrariness of the gravitational field. 

Hence let us now assume that this gravitational degree of arbitrariness is removed by 
specifying completely the gravitational field. Requiring that the space is Galilean has 
this effect, as well as using Newton's equation together with the appropriate boundary 
conditions. 

This done, all invariants involving Tij and its derivatives are removed from our system. 
In order to satisfy condition 01, we need to take the involutive seeds as 

Mij-fi, Ki-fl 

with the as the smallest index. The condition R is satisfied for a general system: it 
suffices to use Mij to solve for ujij, Ki to solve for 9i, and Ki-^o to solve for r. Hence, the 
degree of arbitrariness is 

n(n — 1) 
2 
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occurring at dimension 1. For the usual dimension n = 3+1, this is 6. The dimension of the 
degree of arbitrariness, together with the form of the involutive seeds we use, means that 
we can specify arbitrarily the time derivative of the rotation and acceleration of the rigid 
body at any given time. The directions of evolution are all of the spatial directions: thus 
the constraints we have added do not have interpretations of proper equations of motion. 

Remark 15. With this framework, deriving the degree of freedom of classical gravita- 
tional field is trivial. Indeed, Poissons' equation is simply 

ra-l 

Tii = source terms, 

i=0 

and consequently the contributions to the degrees of arbitrariness is 2, giving the degree of 
freedom of the classical gravitational field 1, similar to the scalar field case. ^ 

3.5 Relativistic rigid flow 

It is perhaps an overkill to use our algorithm to study Newtonian rigid flow. But when it 
comes to relativity, the situation is vastly different. Shortly after the birth of the special 
theory of relativity, Born [3j defined the concept of rigid flow, as the light-like flow that 
drags the orthogonal projection of the spacetime metric onto itself. In other words, for 
a rigid flow, the infinitesimal distance between neighbouring pairs of particles measured 
orthogonal to the worldline of either remains constant. Such a definition remains valid 
when we are in a general Riemannian spacetime. 

Whereas rigid flow in classical physics are simple, thanks to the absolute splitting of 
space and time, in relativity the situation is very different due to the close interplay between 
space and time. Thus, there is no intuitively obvious answer to the following two questions: 

1° What is the degree of arbitrariness, including the degree of arbitrariness of the 
underlying space, of a relativistic rigid flow? 

2° What is the degree of arbitrariness of a relativistic rigid flow in a specified spacetime? 

In the moving frame formulation, a rigid flow in a Riemannian spacetime is completely 
specified by the moving frame 

duJi = —TTij A LUj, 

dujQ = —KiUQ AuJi — MijUJi A loj, 

diTij = —Tlik A TTkj + \Sijkll^k A UJU 

where ujq is aligned with the flowline, oji the orthogonal directions, and the connection for 
the Riemannian spacetime is given by 

(3.10) UjQi = -LOiO = Kil^Q- MijUj, OJij = TTij + MijLVQ. 

The derivation is slightly non-trivial: see [21j. As in the Newtonian case, Ki and Mij has 
the physical interpretations of acceleration and rotation (usually called vorticity in fluid) 
respectively. Sijki is the Riemann tensor for the "reduced space" , analogous to the Riemann 
tensor fixed to the moving body in classical rigid motion. 
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The fundamental invariants are thus Ki, Mij and Sij^i- Sijki has the same defining and 
Bianchi relations as the usual Riemann tensor. Mij has defining relation 

Mij = -Mji 

and we can derive the following additional Bianchi relations 

(3-11) M^ij.j,] = 3M[ijK^], Mij-o = -K[i.j], 5'ijfei;0 = 0. 

Again we have abundance of invariants so that the condition R is satisfied. We can thus 
take the involutive seeds to be 

Mij-k {i > j, i> k) Ki-o Ki-j Sijki-m > J, k> I, i>k, j >l, k> m), 

and the involutive ordering is the ordering such that is the greatest index: this is different 
from the classical case, where we are forced to take the time index to be the smallest. The 
degree of arbitrariness is thus the number of seeds of the form Ki-Q, which means that the 
system has degree of arbitrariness 

n — 1 

occurring at dimension n. 

Remark that we have not yet specified the geometry of the Riemannian space for which 
the flow takes place. In the classical case, we see that the degree of arbitrariness of the 
rigid flow is completely eclipsed by the degree of arbitrariness of the underlying geometry, 
but here, a general Riemannian geometry has degree of arbitrariness n(n — l)/2, which for 
n > 3 is greater than the degree of arbitrariness n — 1 that we have obtained. This means 
that not all Riemannian spaces with n > 3 admit rigid flows (Or, more precisely, almost 
all Riemannian spaces with n > 3 do not admit rigid flows. For n = 2, the situation is 
different, and there always exists rigid flows for any Riemannian space: see [21j). 

What are these n — 1 degrees of arbitrariness? We see that we can freely specify the 
time derivative Ki-Q everywhere. By taking Ki as known functions and use our algorithm 
again, we can show that it is also consistent to freely specify K^. However, it is in general 
not consistent to specify the geometry Sijki of the reduced space independently in space: 
in relativity, we cannot start with a "shape" that we want to move around in spacetime 
and drag it along to obtain a rigid flow, in constrast with the Newtonian case. 

We now come to the second question, namely the degree of arbitrariness in a specified 
spacetime. Using ()3.10p . the invariants Rfj,upX of the total space is related to our invariants 
as follows 



(3.12) 



Rijkl = Sijki + MiiMjk - MikMji - 2MijMki, 
RijkO = Mij.k - MjkKi + M.kKj + M^jKk, 
RoiOj = MikMjk — K(i.j) — KiKj. 



To completely specify a spacetime means simply to take R^ypx, as well as all of their 
derivatives, as known functions. The first relation can be used to express Sijki completely 
in terms of Mij and Rijkh the second can be used to express M^j-k completely in terms of 
RijkO, Mij and Ki (the reader can check that this is really the case, by verifying that the 
symmetries of the equation does not force us to miss any terms), and the third equation 
expresses completely in terms of Roioj, Ki and Kj. By considering Riju-o, we see 
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that in the general case Mij-Q is completely specified in terms of Rijki-o- in more details, 
the first equation implies 

where a bar over i means that there is no summation between i and i, but otherwise i = i. 
This implies 

and if Mjj 7^ for all which is what we mean by the general case, Mij-Q is completely 
determined by R^j-jfl^ which by our assumption is a known function. Now we see that Sijj^i, 
Mij-k, Mij-Q, K(^i.j^ are all completely determined. By (|3.1ip . is as well. Differentiating 
the second equation of (|3.12p and taking those linear in ujq shows that in the general case, 
Ki-Q is completely determined. By our definition in section 12.41 the invariants are covered 
by the empty set, and hence this system has no degree of arbitrariness at all. This means 
that the system either has solutions depending on constants, or has no solutions at all. 
This is in accordance with our remark earlier that, since the degree of arbitrariness of 
the system without specifying the geometry is less than the degree of arbitrariness of a 
general Riemannian space, not all configurations admit solutions. In fact, if the space is 
homogeneous, i.e., R^upx is a single constant, then a generalisation of the Herglotz-Noether 
theorem applies, and all solutions are rotating Killing vectors in the space. 

What about the non-general case? For simplicity, we study only the case where the 
total space is homogeneous. We need to consider only the case where Mij = for all 
z, j, since otherwise by using the right action of the principal bundle we can transform the 
system locally into an equivalent one for which Mij 7^ for all i,j. Then, the equations 
(IXT^ become 

Rijkl Sijki, 
RijkO = 0, 

RoiOj = —K(i-j) — KiKj. 

Immediately we see from the second equation above that, for homogeneous space, unless the 
space has vanishing curvature, i.e., Minkowski space, there is no singular solution. Hence 
assume that we are in Minkowski space. Then 

Sijkl = 0, = -KiKj, 

which together with = 0, forces us to take the involutive seeds as 

Ki;Q. 

Let us now come to the verification of condition R. In the general case, we use Ki to solve for 
uji, Ki-Q to solve for ujq. As for ojij, we use the series of invariants iCj;oo, -f^j;0005 -^i;00005 ■ ■ ■ ■ 
this is possible, since all of these transform under SO{n — 1) and are independent. 

The degree of arbitrariness is thus n — 1, occurring at dimension 1. Intuitively, this is the 
motion of hyperplanes in Minkowski spacetime, the hyperplane always being orthogonal to 
the worldline of any one of its point. The n — 1 degree of arbitrariness is the acceleration 
of any one of its point as a function of time. 



Historical Note. The reference [H] has a good summary of the history and proof of the 
Herglotz-Noether theorem mentioned above, except that it missed out the contributions 
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of Estabrook and Wahlquist |16 t l28 t [25]. which generahsed the theorem to conformany 
flat spacetime. The Herglotz-Noether theorem and its generaUsation by Estabrook and 
Wahlquist are both resticted to 3 + 1 dimensions since the problem was phrased in a form 
that is unsuitable for calculations in arbitrary dimensions: this is especially the case for 
the generalisation, since the dyadic method used by Estabrook and Wahlquist do not apply 
in higher dimensions at all. The present author, on the other hand, has generalised them 
further to all dimensions and all conformally flat spaces, by formulating Born-rigid flows 
as structure-preserving submersions, of which the present example is a very short preview. 

In fact, structure-preserving submersions, when the "structure" in question is the 
Riemannian metric, as it is the case here, is just the familiar problems of Riemannian 
submersions (the reference [1] has a derivation of the structural equations, see also the 
references contained therein for further references on Riemannian submersions). But we 
can also have more exotic structures that are preserved: for example, if a Weyl structure is 
preserved in a Riemannian spacetime, in the codimension 1 case the equations describe pre- 
cisely flows with vanishing shear but possibly non-zero expansion [21]. Needless to say, the 
present algorithm is very useful in such studies since in all cases the degree of arbitrariness 
can be very easily obtained. ^ 

4 Proof of the algorithm 

Now we prove our algorithm. The proof consists of three parts: 

1° Construct an exterior differential system for all the invariants of the system that 
must be satisfied by any solution of the system; 

2° Calculate the degree of arbitrariness (the number of free functions) of the general 
solution to this exterior differential system; 

3° Prove that the solution that we obtained is compatible with the original exterior 
differential system, i.e., the structural equations of the moving frame, hence the 
degree of arbitrariness we obtained really is the degree of arbitrariness of the system 
we are interested in. 

But before that, we need to state the mathematical theorems that we will use. It turns 
out that, for stating these theorems, we need to define a long train of concepts that are 
not familiar to the usual physicist. We put these constructions and statements of theorems 
in appendix [XI and from now on we will use the constructions and results without further 
explanation. See [ll[9l[22] for more in-depth exposition on the theory that we use. 

4.1 The exterior differential system for the invariants 

First we shall assume that the condition R is satisfied, and there are no additional functions 
added to the system. We now consider all the algebraically independent invariants to be 
independent variables. A subspace formed by a subset of these invariants will be the space 
on which we work. 

Since two functions /, g, are functionally independent if and only if their differentials df, 
dg are linearly independent, for systems satisfying the condition R, we can find a full set of 
independent invariants such that the co-frame uj^ is solvable in terms of functions of these 
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invariants. Thus, our original structural equations are now considered differential systems 
whose variables are the invariants: we do not need any coordinates. (More precisely, in 
Cartan's theory of equivalence [8l[18l[22l[23], it is the functional inter-dependence of the 
invariants that specify completely the solution of the system, not the explict analytic form 
taken by these invariants.) 

Now for any invariant Ii-j, we have the expansion 

which by our assumption, only Ij-jk may contain new algebraically independent invariants. 
We will take such equations for a certain set of invariants to be the differential system for 
the invariants, taking care to constrain them by all the Bianchi relations of the systems. 
The forms w^, oja are now considered to be nothing more than shorthands for some linear 
combinations of the differentials of the invariants. This is an exterior differential system 
with independence condition: the independent one- forms are exactly the one- forms Wj, 
cJq, which gives an implicit independence condition for the invariants themselves. For the 
moment we will ignore the original structural equations. 

Our differential system is closed by adding the equations obtained by exterior differen- 
tiation: 

(4.1) = (fli-j = d{Ij.jkUJk + Ci-jaUJa)- 

Observe that, if the invariants Ij-jki for all I are included in our sets of invariants, then 
such equations are identically satisfied: parts of the algebraic relations satisfied by Ii-jki 
are exactly those that ensures this equation holds. Thus, for the two- form equations in 
our differential system, we only need to consider those coming from the differentiation of 
one-form equations consisting of the most number of derivation indices. 

Expanding ()4.ip and using the equations of the differential systems themselves, we get 

(4.2) = dli-jk A Wfc + • • • 

where the dots denote two- forms formed with the independent one- forms. By the previous 
remark, such equations are non-vacuous only in the cases where derived invariants of Ij-jk 
are not included in our system, in which case dlj-jk are now independent one-forms, and 
hence they are the only two-form equations in our differential system. 

Observe that our definition for a set of involutive seeds is tailored such that, for the 
involutive seeds and all the invariants covered by them, a set of independent and spanning 
one-forms dij-jk can be taken exactly as those of the derivatives of the involutive seeds: the 
condition II ensures that the system includes all the invariants that must be included, and 
12 ensures independence of the seeds. The case where the set of involutive seeds is empty 
occurs if and only if when including enough differential invariants into our system, the sys- 
tem is Frobenius-integrable (or inconsistent), and hence contains no degree of arbitrariness 
by Frobenius theorem. 

4.2 The degree of arbitrariness of the system of invariants 

Now we apply the Cartan-Kahler theory to our differential system. We will use version 2 
of the involutivity theorem, and hence what we are interested in is the reduced characters 
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without taking into account of the maximal rank condition. As the involutivity test does 
not need s'q, and all the Cartan characters except s'q depend only on equations formed with 
two-forms and higher order forms, we only need to focus on (|4.2p . The reduced characters 
in this case are calculated by feeding the two-forms an independent vector direction so that 
it becomes an one-form equation, and ignore all the independent forms. Thus, due to the 
form of (|4.2p . we see that, for any ordering of the one- forms ujk, which without loss of 
generality we take to be 1, 2, ... , n, and label the other one-forms Ua, a = n + 1, . . . ,m, we 
have 

s'l is the number of independent invariants 

S2 is the number of independent invariants Ii;j2] 

is the number of independent invariants Ii-^jn] 
s'^ for n < k <m\s equal to 0. 

In the above, when we say for example "the number of independent forms Ii-jk", we mean 
that the Ij-jk are not only independent among themselves, but also with respect to all Ij-ji 
for I < k. Thus, for any system of involutive seeds satisfying in addition condition 01, in 
particular, for any system of involutive seeds together with an involutive ordering, s- is just 
the number of involutive seeds with the last index equal to i. 

However, the second version of involutive test concerns also the pseudo-character. We 
now show that the pseudo-character either vanishes or is equal to the non- vanishing s'^. 
Observe that, for every independent invariant that is not a seed there is an independent 
one-form equation. Hence s'q is the number of independent non-seed invariants minus m, 
as m of these have been used implicitly to solve for the w^. The above counting, on the 
other hand, shows that the number of seeds is given by X^j^^ m^i^ 

s'o + s'i-\ \- s'jn = M - m, 

where M is the dimension of the space formed with all the independent invariants included 
in our system, showing that the pseudo-character and the last character are equal. 
Using the second version of the involutivity test, if we can show 

(4.3) s'l + 2s^ + ■ ■ ■ + ns'^ = N 

where N is the number of free parameters for an integral element, then the system is 
involutive. To calculate the number of free parameters, we expand dij-jk, which hitherto 
has been considered to be independent, in terms of cok, ^a- But this expansion is just the 
definition of the invariants Ii-ju- 

dlljk = Il;Jkl^l + • • • 

where the omitted terms contain no free parameters, by our assumption. Hence N is the 
total number of algebraically independent first derivations of the involutive seeds. Since we 
have included enough invariants in our original set, we see that the Ii-jki are only subject 
to the derived relations and the generic relations. Hence an independent and spanning set 
of Ii-jki can be taken as those Ii-jki for which k > I (using the generic relations) and Ii-jk 
is an involutive seed (using the derived relations). 
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Now, for a system of involutive seeds together with an involutive ordering, this, together 
with the condition 02, means that 

The number of indep. Ii-jki is the number of seeds Ii-ji, Ii;j2, ■ ■ ■ , Ir,jn', 
The number of indep. Ii-jk2 is the number of seeds Ir,j2, ■ ■ ■ , Ir,j,n', 

The number of indep. Ii-jkn is the number of seeds Ir,j,n- 

Thus, (j4.3p is satisfied, and the system is involutive. Our interest hes in the degree of 
arbitrariness, i.e., the last non-zero Cartan character, and this is just the number of seeds 
with largest the last index. If this index is k, then the degree of arbitrariness = has 
the interpretation that the general solution of the system of invariants depends on Sk free 
functions of k variables. 

A last very important thing that we need to note is that, by applying the theorem on the 
complete prolongation of Pfaffian systems, the system of invariants, which we have found 
to be involutive, can be prolonged at will, with the degree of arbitrariness unchanged: in 
particular, our system can be taken to include any given invariants. 

Note that since the coordinates we use are provided by the invariants themselves, we do 
not have to worry about any degree of arbitrariness coming from the "freedom" in choosing 
the coordinates. This is what is meant by saying that our system is "coordinate- free" : the 
price to pay being the condition R. 

The claims of the "directions of evolution" follows directly from the interpretation of 
the Cartan-Kahler theorem. 

4.3 The original structural equations 

We have thus obtained a solution to the differential system for the invariants. With this 
solution, the one-forms and uja can be expressed in terms of any set of invariants which 
we take to be the coordinates. It remains to show that the one-forms ujk and Ua thus 
constructed satisfy the original structural equations of the problem, which we have ignored 
so far. Since the differential system we use to obtain the solution respects the independence 
condition, we can find a set of functionally independent invariants whose number is equal 
to the number of independent forms. For I^, our differential system include the equation 

(4.4) dl^, = I^,,uUJu 

where our notation has changed so that uj^ includes both cok and coa- In this notation, our 
original structural equations are written 

which are not included in our differential system, and the derivative is written as 

which can all be included in our differential system due to the possibility of unlimited 
prolongation, and the invariants are chosen among the c^^^a and its derivatives. 
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As both dl^ and oj^ are co-frames for the solution integral variety, the matrix I^^^, is 
non-singular. Deriving ()4.4p . we get 

(4.5) = dlf,^i, Au},^+ If^^udoJu = I^i,u\^^v A lox 

But as now we can take any invariant in the system we have used to solve for the degree 
of arbitrariness, is related to for example, as 

we have 

= d^c^yx = c^i/A,75W7 A a;5 + c^j^yx^-fC^pS^p A ujs, 

where uj^ps are the fundamental invariants. It is essential to note that this equation is a 
consequence of our differential system for the invariants and we do not need to assume the 
original structural equations: indeed, the algebraic relations for c^yx,-iS in our differential 
system are required to include those relations using which this relation holds. Thus, ()4.5p 
can be written 

(4.6) Ip,,v{duJu - Cupx^p A wa) = 0. 
As Ifj,^i, is non-singular, 

du),, - c^pxUJp A = 0, 

and the original structural equations are satisfied by our solution. Thus the proof of our 
algorithm is complete for the case where condition R holds and there are no additional 
functions. 

4.4 The case with non-maximal number of invariants 

The problem with the case where the condition R fails is that, first, it is impossible to solve 
all of in terms of the differentials of the invariants, and thus it is not possible to form a 
differential system for the invariants in the way that we did for the non-singular case, and 
second, even if the first difficulty is somehow overcome, in ()4.6p the matrix I^^j^ is singular, 
and thus there is no way to ensure that the original structural equations are really satisfied. 

Since in this case, among the differential of the invariants, there are only p, < p < m 
linearly independent ones where m is the total dimension, there are also only p functionally 
independent invariants. To make progress, we first consider the equivalence problem for 
such a system. The method of equivalence is briefly summarised in appendix [Bj 

For an equivalence problem formulated on the manifold M of our system, let us consider 
the Euclidean space C of the invariants Cijk up to sufficiently high order so that all function- 
ally independent invariants are guaranteed to be included. The exact order is immaterial, 
that we have considered a lot of redundant variables does not matter either. Then, for every 
concrete system, we have a classifying map T : M — > C. This map defines a classifying 
manifold in the space of invariants. Using the classifying map and manifold, we can state 
the equivalence condition as follows: if two points P and P of two systems defined on M 
and M map to the same point in the classifying manifold (the classifying manifolds of the 
two problems are identified in the obvious manner), then the two systems are equivalent 
at the points P and P; if this condition holds for all points in the open sets S C M and 
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S C M, then the identification of S and S that makes the condition holds provides an 
equivalence of the two systems. 

For our present purpose, we are interested in the case where M = M. Let us focus on 
the image of the classifying map T : M ^ C. At regular points, this image is a submanifold 
of C: its dimension is exactly p. li p < m, then the pre-image of a point Q G C is non- 
trivial: T^^{Q) is locally a m — p dimensional submanifold in M. But if we set up the 
equivalence problem in trying to deduce the equivalence of M with itself, but identifying 
a point P in M with a nearby point P' in a neighbourhood, we see from the reasoning of 
the previous discussion that as long as P and P' are in the same pre-image, i.e., as long as 
T{P) = T{P'), the equivalence problem has a solution. As P and P' can be connected by a 
path not going out of the pre-image, we see that this self-equivalence is actually a symmetry 
Tinder a fintc dimensional Lie group, and m — p gives the dimension of the symmetry group 
of the problem. (If p = m, then no non-trivial continuous symmetry group can exist.) 
Thus, for systems with a non-maximal number p of invariants, we can find m — p vector 
fields Vq, corresponding to the symmetry group, for which 

C^Un = 0, v{cijk,lm...) = 0. 

Now suppose we take a submanifold N of dimension p transverse to all of the the vectors 
Vq, and we find a solution of our system (i.e., a functional dependence of the forms u;^ and 
the invariants c^i^p on m coordinates satisfying the structural equations) valid in an 
infinitesimal neighbourhood of N, then using the system of vector fields Vq, this solution 
can be extended to the whole space. Also, any solution valid for the whole space M, when 
restricted to such a transverse submanifold N, will also satisfy the structural equations: 
the equations that are to be satisfied now are just the pullbacks of the equations on the 
total space, and exterior differentiation commutes with pullbacks. On the other hand, 
the converse is in general not true: for simplicity, let xi,X2, . . . ,Xp be the coordinates on 
the submanifold A'', and Xp+i, . . . ,Xm be the transverse coordinates. Then our structural 
equations are written in the differentials of these coordinates. That we have a solution on 
N means that, for the structural equation, 

dujp{dx^,dxp) = CpupUJy Aijjp{dx„,dx0), {a, f3 < p) 

is satisfied on at points on A^. We can even show that for certain coordinates and for 
a,l3 > p, this equation also holds, due to the action of the Lie group. However, in general, 
when a < p and P > p, there is no reason that this equation will hold. Hence, that the 
structural equations are satisfied when pulled back onto any submanifold A^ thus chosen is 
a necessary, but in general not sufficient condition for a solution of the original equations. 

Now let us return to our differential invariants and involutive seeds. First observe that, 
if Ii;Jk is an involutive seed, then cok depends on the differentials of the invariants: indeed, 

dli-j = Ii-jkUJk 

and due to the requirements of the involutive seeds, when we prolong the problem, 

dli;jk = Ii-jkk'^k "I 

where Ij-jkh is an independent invariant. This shows that for all solutions, we can find a 
submanifold A" such that all the forms lj^ occurring with a differential of the involutive seed 
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in equation (|4.2|) remain independent. On other other hand, clearly if for the system (|4.2|) . 
if we take any m — p forms other than those occurring explicitly with the involutive 
seeds in ()4.2p to be forms written in terms of the invariants, with arbitrary functional 
dependence, we can now solve the remaining p one- forms o;^, which contains the cok, in 
terms of the invariants. Due to the form of the equation (j4.2p . for any such choice, the 
Cartan characters for this system, with p independence conditions, are the same, and by a 
reasoning exactly the same as in the maximal rank case we see that this system is involutive, 
with degree of arbitrariness given by the last non-zero Cartan character. Observe also that 
if for a system of involutive seeds and ordering, if the maximal last index is k, then p > k: 
this can be seen easily from (j4.2|) as well. This means that the degree of arbitrariness 
always occurs at a dimension < p. 

Using the Lie group action, such a system of values of the invariants can be extended 
to the whole space. But as remarked earlier, there is no guarantee that after extension all 
of the structural equations will be satisfied, hence in this case we have obtained only an 
upper bound. 

Finally, let us remark that for the following special case, the upper bound is realised: 
the structural equations reads 



where Cq^^ are constants, and all of coi can be solved in terms of Cijk and their invariants: 
it suffices to first ignore the first set of equations and obtain the degree of arbitrariness 
for this system. This degree of arbitrariness is realised since the first set of equations is 
consistent (otherwise the Bianchi relations will have something that reduces to 1 = 0), and 
from the theory of Lie groups we know that there exists a Lie group satisfying the first set 
of equations. The solution space is then the product space formed by the Lie group and the 
solution whose degree of arbitrariness we know. As an easy consequence, for the following 
system the upper bound is also realised: 



where except for the modulus part, the assumptions are the same as before. 

4.5 Additional functions 

We only consider the case satisfying condition R, as our algorithm states. In this case, 
we simply adjoin the additional functions and their derivatives to the system (j4.2p . and 
the result is obvious. The condition 13 is necessary since only those functions satisfying 13 
will appear in the written out part in the right hand side of (j4.2p . for which the Cartan 
characters are calculated. 
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A Cartan-Kahler theory, involutivity and prolongation 



Exterior differential systems. Let a exterior differential system be defined on a manifold 
M, formed by 

functions: fi, ■ ■ 

one-forms: ^i,-- 

two-forms: , 

three-forms: ©i^^ 

We also assume that that the given system has already been transformed into a closed 
system: 

d9i = (mod O'^f), etc. 

so all relations that can be obtained by differentiation have already been incorporated. A 
solution (integral variety) of this exterior differential system is a submanifold of M such 
that all the forms in the system vanish when restricted to this submanifold. 

Characteristics. Before talking about Cartan-Kahler theorem, we need to talk about 
characteristic directions. A direction given by the vector v is characteristic if 

V J G = for all G in the differential system, 

i.e., characteristics directions are the ones having the property that if an clement contains 
this direction, then it is automatically an integral element. To obtain the integral vari- 
eties of the characteristics, define the characteristic system C associated with our original 
differential ideal, consisting of all one-forms i? such that 

??(v) = for all V G V. 

Of course, not all systems contain characteristic directions. 

The characteristic system is important when we are solving the Cauchy problem. Sup- 
pose we are integrating a differential system and we want to go from dimension k to 
dimension k + 1. We need to specify the initial data on the integral variety of dimension 
k which we have already found. Biit if this variety contains a characteristic variety, then 
we know two things: first, the data along these varieties are well-defined when we know 
them at any point on them, hence there are consistency issues when specifying the initial 
data; second, we need specify extra functions in order to effect the integration. Hence for 
the Cauchy problem, specifying initial data on varieties containing characteristic variaties 
is problematic and in general inconsistent. 

In the following, we assume that any vectors and any submanifold that we talk about 
do not contain characteristic directions. 

Integral elements; Cartan characters. A first step in constructing integral variety is 
constructing integral elements: these are subspaces of the tangent space at a single point 
which satisfies the equations. We construct the integral elements recursively. Obviously, it 
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suffices to restrict our attention to integral elements over the submanifold defined by the 
vanishing of the functions in our differential system. 

First let us consider integral elements at a single point. Every one-dimensional integral 
element v must satisfy 

^^(v) = 0. 

Hence the space of integral element of dimension 1 is given by a linear equation: Aij are 
constants at the points we consider and any solution Vj forms just a vector. Now suppose 
we already have a determined 1 dimensional linear element Vi and we would like to extend 
it into a 2 dimensional linear element: this amounts to finding another direction V2 = Vjlj. 
Wc need to ensure AijVj = 0: this direction must itself be a solution to the one-dimensional 
problem. But there is more: suppose the 2-forms in the system are 

then we must have 

AjkVjVk = 0. 

Vj is data already given to us: this, together with AijVj = 0, forms a linear system whose 
solution space gives all possible directions extending vi. In general, given a set vi, . . . , Vp 
forming an integral clement of dimension p, the extension to dimension p + 1 is obtained 
by a suitable linear system on the free tangential directions. This linear system is deduced 
from all the exterior forms in the system of dimension up to p + 1. As we go up in the 
dimension of the integral element, the dimension of the solution space will decrease: for 
every step it will decrease by at least one: for example, at the second step, vi and V2 must 
be independent for them to constitute an extension of vi . Thus we will eventually come to 
a dimension where the integral element can no longer be extended. 

The integral varieties we are looking for have something to do with the ranks of the 
various linear systems we have just described. But we cannot ensure that, for all choices of 
the point P in M, the rank of the system AijVj remains the same. Thus let us call a point 
P an integral point if it is on the algebraic variety defined by the vanishing of functions 
/i = in the differential system. An integral point P is generic if, at the point, dfi = fijdx^ 
has maximal rank in a neighbourhood. This ensures that we are not dealing with some 
pathological algebraic variety. 

On a generic point, let us find the one-dimensional integral elements v. If the rank of 
Aij defining the integral element is maximal in a neighbourhood, then the integral element 
is said to be ordinary and this generic point (not the integral element!) is said to be a 
regular point. The rank of the matrix A^j is called the zeroth Cartan character (or simply 
zeroth character) of the system at the point P and is denoted Sq. 

Now given an ordinary one-dimensional integral element Vi at a point, let us try to 
extend it by one dimension. This element is said to be regular if the rank of the system 
AijVj, AijkVjVk, where Vj is now given, is maximal in a neighbourhood. The solutions for 
such a maximal-rank system vi, V2, are said to be ordinary integral elements. Since the 
rank of the system AijVj, Aij^VjVk cannot be less than the rank of a part of itself, AijVj, 
we denote the rank as sq + si. The integer si is called the first Cartan character. 

The general pattern should be clear: for the linear system defined at a particular di- 
mension p, the maximality of its rank defines the ordinary elements at this dimension, and 
the regular elements at one less dimension. Then, as we have remarked that the dimension 
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of the integral elements cannot increase indefinitely, we will come to a dimension p where 
we have ordinary elements that cannot be further extended, hence there are no regular 
elements at this dimension. The rank of this final system is 

So + si H H Sp, 

and for such a system we have a set of p + 1 Cart an characters. 
After this long train of definitions, we can finally state 

Cartan-Kahler Theorem. On a manifold of total dimension n, for a given integral variety 
admitting ordinary integral element of dimension k at a point, we have 

(A.l) so + siH \- Sk<n — k, 

equality holds above if and only if k is the largest such dimension. If all data used to define 
the exterior differential system are analytic functions, then the above condition is also 
sufficient: i.e. there will exist a regular k-dimensional integral variety for the differential 
system. 

Let us also define the Cartan pseudo-character by the formula 

So + si H \- Sfe-i + cTk = n - k. 

It is always non-negative by the theorem, and if k is the maximal dimension, cjfc = s^. 

Remark 16. The sufficiency part of the Cartan-Kahler theorem requires analyticity. For 
physical applications, should this worry us? 

First, let us note that it is sometimes possible to substitute the Cauchy-Kowalewski 
theorem for some stronger existence theorem in the theory of partial differential equations: 
thus the Cartan-Kahler theorem actually holds in the smooth setting (as opposed to the 
analytic setting) under some circumstances. In particular this is true for involutive hyper- 
bolic systems, quasilinear systems, and some Pfaffian systems. Most of the equations that 
arise in physics fall into one of these classes, so we actually have guaranteed C°° results. 

Second, as the necessity part does not require analyticity (neither does the "degenerate" 
case where the theorem reduces to Frobenius theorem), we will never miss any solutions 
by applying the Cartan-Kahler theorem in cases where the C°° version does not hold: the 
worst we do is claiming the existence of solutions where there are none. But as the space 
of analytic functions is dense in the space of continuous functions, there is no physical 
ground for claiming that a function is differentiable to a high order but is not analytic: 
physically we cannot really distinguish, say, data from C"^ data since the error in our 
knowledge would already make them indistinguishable. Furthermore, if we are really in a 
situation where such issues becomes significant, for example due to the discontinuity of the 
underlying space, maybe with origin in some kinds of "quantum" effects, there should be 
a lot more other things to worry about than analyticity. ^ 

Interpretation of the characters. For our purpose, the following interpretation of the 
theorem is important: the regular fc-dimensional integral variety guaranteed by the theorem 
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can be integrated by specifying 

si + S2 + • • • + Sfc-i + CTfe arbitrary functions of x^, 
S2 + ■ ■ ■ + + CTfe arbitrary functions of x^, x^, 

Sk-i + CTfe arbitrary functions of 

(Tfe arbitrary functions oi x^,x^, . . . , x'^~^,x^ . 

where X ^ X 5 . . . IS cL coordinate system system for the manifold. Furthermore, the last 
non-zero (pseudo)-charactcr is not affected by our choice of order of integration, the form 
we used to write down the differential system, etc. (the characters before that are affected 
by such choices). Obviously, for a differential system, the last non-zero (pseudo)-character 
is its degree of arbitrariness when we are interested in forming fc-dimensional integral vari- 
eties, and the index of this (pseudo)-character is the dimension at which this degree of 
arbitrariness occurs. 

Independence conditions; involutive systems. The solutions guaranteed by the Cartan- 
Kahler theorem has a problem: it does not care for our requirement of independent vari- 
ables. To make progress, let us define a differential system with independence condition as 
a differential system for which we require that solutions must keep certain one- forms wi, 
U2-, • • • , ojm independent. This amounts to 

uji A LUj /\ ■ ■ ■ A Uk 0, ... ,k all distinct 

on solutions for any choice of any numbers of coi, coj, . . . . For the coordinates xi, X2, • • • , 
Xm, we just take the forms to be dxi, dx2, ■ ■ ■ , dxm, and the vectors spanning the integral 
elements we look for must be of the form 

(A.2) ^+ ± „ = 1,....™. 

" i=m+l * 

Immediately, we see that we must not end up with equations of the form 
(A. 3) = = Aij^^^k oJi Aojj A ■ ■ ■ Auk- 

Unless Aij_,,k = 0, no solution will satisfy the independence conditions. Such terms are 
called essential torsion. To proceed in such cases, we need to add to our differential system 
the algebraic equations 

■^ij...k — 0- 

Such equations might have no solution: for example if A^j j. are non-zero constants. We 
say that such a differential system is incompatible. 

Now consider only systems with no essential torsion. Our aim is to obtain solutions 
satisfying the independence conditions as general solutions of the problem without inde- 
pendence conditions. Hence, we will call differential systems for which the constraints on 
ordinary m dimensional integral elements do not require any linear relations among the 
independence conditions wi , . . . , ujm involutive systems. This does not mean that all or- 
dinary m dimensional integral elements satisfy the independence conditions: it means that 



37 



almost all satisfy, and almost all integral elements spanned by the vectors of the form (jA.2p 
are ordinary, since the conditions for otherwise are both equality conditions. For involut- 
ive systems, we can obtain the integral varieties we want by applying the Cartan-Kahler 
theorem to ordinary integral elements of the form ()A.2p . 

Reduced characters; Cartan's tests. The requirement that none of the constraints on 
integral elements can involve any relations for Wj means that, in calculating the Cartan 
characters, we can ignore all directions that correspond to Ui For example, we have two 
forms 

17" = A^oji A ooj + B^oji A vr, + CS,vr„ A vr, 

and in calculating the Cartan character, the first direction is chosen as the vector Ii dual 
to LOi. With this direction, we have 

The requirement that we must not have any constraints for the independent directions 
(i.e., no relations of the form CjWj = 0) means that the rank of the system {Afj, is the 
same as the rank of the system {B"^). Notice we do this only after using the already found 
tangent direction. 

Let us call the numbers 

*0) '^2) ••• *m-l 

which are calculated by omitting all terms corresponding to cjj in the calculation for Cartan 
characters the reduced Cartan characters: this definition holds for both involutive and non- 
involutive system. The reduced pseudo-character cr^ is defined as 

s'o + s'l + • • • + s^-i +<7'^ = n-m. 

Involutive system has the normal characters and reduced characters equal. Conversely, 
if all reduced characters are equal to the non-reduced counterparts, this implies there is no 
relation among the independent variables, and hence the system is involutive. We therefore 
have a necessary and sufficient condition for involutive systems: the equality of the reduced 
and normal characters. 

It is troublesome to calculate two sets of characters, especially the non-reduced ones. 
The following theorem gives a simple way to tell if a differential system with independence 
condition is involutive: 

Cartan's Involutivity Test. A differential system with independence condition is involutive 
if and only if 

(A.4) N = s[ + 2s'2 + --- + {k- l)4_i + kal 

where N is the total number of free parameters occurring for the ordinary k-dimensional 
integral element. 

Calculating reduced characters is easier than calculating the normal characters, but we 
still need to ensure that the integral element we specify is generic, meaning that the rank 
of the various systems are maximal. Now assume that we calculate the reduced characters 
by ignoring the step of checking the rank condition. Then we still have 
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Cartan's Involutivity Test, Version 2. A differential system with independence condition is 

involutive if 

(A.5) iV = s'l + 24 + • • • + (fe - l)4_i + kal 

where there various reduced characters are calculated without checking the rank condition, 
and N is the total number of free parameters occurring for the ordinary k- dimensional 
integral element. 

Note that there is no "only if" for this easier test. In practice, for calculating the 
number of free parameters, we take the differential one-forms Wa that are not required to 
be independent in our system, and write 

where uji are the independent forms. Then substitution into the system gives algebraic 
constraints on A^j. The components that are not constrained gives the number of free 
parameters. 

Prolongation. Actually, when formulating the original problem, we can also take the 

some of the parameters A^i occurring above as variables. Doing this is called effecting 
a prolongation of a differential system (by addition of new variables). If we take all of 
the parameters as new variables, we are effecting a complete prolongation. The Cartan- 
Kuranishi theorem states that for well-behaving systems (which we will not make precise 
here), after a sufficient number of complete prolongations, we either arrive at an involutive 
system or an incompatible system. For us, on the other hand, we need the following 
theorem: 

Theorem on the complete prolongation of Pfaffian systems. If a differential system consists 
only of functions, one-forms, and two-forms that are generated by differentiating the one- 
forms (i.e., a Pfaffian system), and the system is involutive, then its complete prolongation 
is also involutive, with the Cartan characters of the new system related to the old system 
by 

sl = Si+ S2^ h S,n-1 + Sm, 

S2= S2-\ h Sm-l + Sm, 

Sm—1 "I" 

Sm- 

This theorem also explicitly verifies the general result that the last non-zero Cartan 
character is an intrinsic property of the system, i.e., not affected by the way of presentation. 

B The method of equivalence 

Our task is to understand, given two co-frames 9i, . . . , On and 61, . . . , On defined on two 
differential manifolds M and M of dimension n, when they are "locally the same" , meaning 
there is a map f : M ^ M such that 

f% = Oi. 



Sm—l 
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We will omit pullback signs from now on and just write 9i = 9i. 

We can treat this equivalence problem as the integrability of the form 

in the space M x M. We are interested in the m-dimensional integral manifolds which are 
transverse, meaning that either 9i or 9i can be taken as the independent conditions. 
As we are dealing with coframes, we can calculate their exterior derivatives 

dOi = Cijk9j A 9k, dOi = Cijk9j A 9^, 

so 

d-&i = CijkOj Ai}k + Cijk"&j /\9k + {cijk - Cijk)9j A 6k 
The existence of integral manifolds then requires 

(B.l) Cjj'fc C-ijki 

i.e., the system must be Probenius integrable (integral manifolds from the the Cartan- 
Kahler theorem will never be of the required dimension, since we already have sq = 2m — 
m = m, and any non-zero Si for i ^ will lower the dimension of the integral manifold). 

The simplest case is where Cijk and Cijk are both constants. If they are equal, then the 
two systems are equivalent, otherwise (including the case where one set are constants while 
the other are not) they are not. 

If Cijk and Cijk are not constants, then we require f*Cijk = Cijk, but for functions under 
diffeomorphism we cannot directly compare them. Worse, even though we can implicitly 
define the submanifold of M x M by the relation Cijk = Cijk, it is far from certain that 
restricted to this submanifold what will happen to 9i and 6i, which we require to be inde- 
pendent. 

To make progress, we treat (jB.ip as a new condition and adjoin it to our conditions for 
equivalence (prolongation). The closure of our system now includes the condition 

dcijk — dcijkj or Cijk-i — c-ijk-^u 
where the coframe derivative for a function h is defined, as usual, 

dh = h-i9i. 

Let us first investigate the case where Cijk and Cijk both contain m independent functions 
among them each. Let us denote these by 

-^1, h, ■■■ Im', h, h, ••• Im, 

and the rest of the quantities are expressible as functions of them. Being functional in- 
dependent, their differentials are linearly independent, or the matrix li-j in the following 
expression is invertible: 

dli = Ii.j9K 

In this case, Cijk = Cijk implies li = li, which in turn implies dli = dli and hence 9i = 6i. 
Observe that we do not even need to check the equality for the terms other than /j. Indeed, 
let us diff'erentiate the above equation. We get: 

= li.jke'' A 9^ + U,jCjkiO'' A 9' 
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so Jj = li implies 6i = 6i and in turn these two equations together imphes the equahty of 
everything else. Hence, in this case, we only need to check 

In the general case, what we do is we use the coframe derivative to differentiate the fun- 
damental invariants Cijk until we get no more functional independent quantities. For 
example, in this way we may obtain a system 

~ ^i-i Cjjk ~ ^ijk-i ('ijkil ~ ^ijk;li ^ijk;lm ~ ^ijk;lmi ^ijk;lmn ~ ^ijk;lmni 

where Cijk-iran introduce no new functionally independent quantities. These conditions are 
obviously necessary. They are also sufficient for equivalent, as it can be easily checked that 
they imply the Probenius integrability of the system 

(^i Gii C-ijk C-ijki ^ijk;l C.ijk;h ^ijk;lm C.ijk;lrm 

and the integral manifold really is transverse. 
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